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Abstract. We present preliminary results on the metastable behavior of a nonequilibrium ferro- 
magnetic system. The metastable state mean lifetime is a non-monotonous function of tempera- 
ture; it shows a maximum at certain non-zero temperature which depends on the strengh of the 
nonequilibrium perturbation. This is in contrast with the equilibrium case in which lifetime in- 
creases monotonously as the temperature is decreasesed. We also report on avalanches during the 
decay from the metastable state. Assuming both free boundaries and nonequilibrium impurities, the 
avalanches exhibit power-law size and lifetime distributions. Such scale free behavior is very sen- 
sible. The chances are that our observations may be observable in real (i.e. impure) ferromagnetic 
nanoparticles. 



Metastability is ubiquitous in Nature, and it often determines the system behav- 
ior. The microscopic understanding of this phenomenon, which is mathematically 
challenging,[l] is therefore of great practical and theoretical interest. In particular, 
metastability is relevant to the behavior of magnetic storage devices. A magnetic 
material may consist of magnetic monodomains. In order to store information on such 
material, one magnetizes each individual domain using a strong magnetic field, defining 
in this way a bit of information. A main concern is to retain the individual orientations of 
the domains for as long as possible in the presence of weak arbitrarily-oriented external 
magnetic fields. The interaction with these external fields often produces metastable 
states in the domains, and the resistance of stored information thus depends on the 
properties of these metastable states, including the details of their decay. On the other 
hand, real magnetic domains are usually impure. The microscopic nature of impurities, 
which shows up in actual specimens as spin, bond and/or lattice disorder and other 
inhomogeneities, quantum tunneling,[2] etc., suggests they might dominate the behavior 
of near— microscopic particles; in fact, they are known to influence even macroscopic 
systems. An interesting issue is therefore understanding the formation of a new phase 
inside a metastable domain which contains impurities. 

Following recent efforts, [3] we study in this paper the simplest possible model of 
this situation, namely a 2d Ising ferromagnet with periodic boundary conditions that we 
endow with a weak dynamic perturbation competing with the thermal spin-flip process, 
which impedes equilibrium. Consider the Hamiltonian function J4?(s) = —J^Uj) s i s j ~ 
hYjjLi Si, where J > (ferromagnetic interactions), Si = ±1 stands for the two possible 
states of the spin at site i of the square lattice, i = l,...,N, and the first sum is over 
any pair of nearest— neighbor sites. The lattice side is L; N = I? . The system 




FIGURE 1. Critical temperature for our nonequilibrium ferromagnetic system as a function of 
p in a first-order mean-field approximation. 



configuration, s = {5,}, is let to evolve in time due to superposition of two canonical 
drives. That is, we chose the transition probability per unit time for a change of ?to be 

co(s^r)=p+(\- P ) p— (i) 

Here s' stands for s after flipping the spin at i, and = J4?(s l ) — J4?(s). The 
Boltzmann constant is set kg = 1 in this paper. One may interpret that rule (1) describes 
a spin— flip mechanism under the action of two competing heat baths: with probability 
2/7, the flip is performed completely at random (as if ?were in contact with a heat bath at 
'infinite' temperature), while the change is performed at temperature T with probability 
1 — p. For p = 0, eq. (1) corresponds to the equilibrium Ising case, which exhibits for 
h = a critical point at T = Tq ~ 2.26917. Otherwise, the conflict in (1) impedes 
canonical equilibrium, and the system evolves with time towards a non— equilibrium 
steady state that essentially differs from the Gibbs state for T . It is assumed that this 
kind of stochastic, non— canonical perturbation for p > may actually occur in Nature 
due to microscopic disorder or impurities, etc. [2] The system shows for p 7^ and h = 
a continuous order-disorder phase transition at a critical temperature T c (p). This critical 
temperature can be calculated in first order mean field approximation as [4] 

TM = zl_ (2) 

1 i»[4+i\/Wi 

Fig. 1 shows T c (p) as a function of p. T c (p = 0) is the Bethe temperature, Tsethe ~ 
2.88547, to be compared with the exact critical value for p = 0, Tq. There is a critical 
value of p, p c , such that for larger values of p there is no ordered phase for any 
temperature. This can be obtained from the condition T c (p c ) = 0, yielding p c = ^ = 
0.15625. In this paper we pay attention to the ordered phase. 




FIGURE 2. Semilog plot of the metastable-state mean lifetime as a function of temperature 
for a system with side L = 53 and, from top to bottom, p = (equilibrium) and p = 0.001 
(nonequilibrium). The data come from an average over 1000 independent runs. Notice the non- 
monotonous behavior of T for the nonequilibrium case. 



In order to characterize the metastable behavior of the model, we measured the mean 
lifetime for values of T and p such that T < T c (p) and a magnetic field h = —0.1. We 
take the initial state to be s,- = +1 for i = 1 , . . . , N. Under the negative field, this ordered 
state is metastable, and it eventually decays to the stable state. In fact, the system rapidly 
evolves from the initial state to a state in the metastable region, with magnetization close 
to +1. After this fast relaxation, the system spends a long time wandering around the 
metastable state, eventually nucleating one or several critical droplets of the stable phase, 
which rapidly grow thus making the system to evolve from the metastable state to the 
stable one. We define the mean lifetime of the metastable state, z(T,p,h), as the mean 
first-passage time (in Monte Carlo steps per spin, MCSS) to m = 0. The simulations 
reported here required in practice using the s-1 Monte Carlo with absorbing Markov 
chains (MCAMC) algorithm, and the slow forcing approximation. [5] Fig. 2 shows 
l{T,p, h) as a function of T for a L = 53 and p = (the equilibrium case) and p = 0.001 . 
Rather amazing, we observe in the figure how the nonequilibrium lifetime exhibits a 
maximum and then decreases as the temperature drops. The maximum occurs at T max (p) 
which depends on the nonequilibrium parameter p. On the contrary, the lifetime in 
the equilibrium case grows exponentially fast as a function of l/T, as predicted by 
nucleation theory [3]. We do not have a simple explanation of this nonequilibrium effect. 
However, it may be emphasized that it has some practical implications. That is, one 
should not blindly decrease temperature when trying to maximize the lifetime of real 
(i.e., impure) metastable magnetic particles, but look for the temperature T max (p) which 
maximizes the lifetime for the typical impurity concentration of the particle. 

The previous observations concern the bulk metastability. However, in real magnets, 
one needs in practice to create and to control fine grains, i.e., magnetic particles with 
borders whose size ranges from mesoscopic to atomic levels, namely, clusters of 10 4 
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FIGURE 3. (a) Time variation of the magnetization showing the decay from a metastable state 
for a R = 30 particle; avalanches are seen here by direct inspection. Time is in Monte Carlo 
Steps per Spin (MCSS), and T ~ 10 30 MCSS. (b) Large avalanche size distribution P(A m ) for 
the circular magnetic nanoparticle and R = 30 (bottom) and R = 60 (top). The second curve has 
been shifted in the vertical direction for visual convenience. 

to 10 2 spins, and even smaller ones. Though experimental techniques are already accu- 
rate for the purpose,[6] the underlying physics is much less understood than for bulk 
properties. In particular, one cannot assume that such particles are neither infinite nor 
pure. That is, they have free borders, which results in a large surface/volume ratio in- 
ducing strong border effects, and impurities. Motivated by the experimental situation, 
we studied a finite, relatively small two— dimensional system subject to open circular 
boundary conditions. The system is defined on a square lattice, where we inscribe a cir- 
cle of radius R; sites outside this circle do not belong to the system and are set = 0. 
We mainly report in the following on a set of fixed values for the model parameters, 
namely, h = —0.1, T = 0.1 ITc and p = 10 6 . This choice is dictated by simplicity and 
also because (after exploring the behavior for other cases) we came to the conclusion 
that this corresponds to an interesting region of the system parameter space, where the 
effects of p and T are comparable and clusters are compact and hence easy to analize. 
We believe that we are describing here typical behavior of our model, and the chances 
are that it can be observed in actual materials. 

The effects of free borders on the metastable- stable transition have already been stud- 
ied for equilibrium systems. [7] In this case, the system evolves to the stable state through 
the heterogeneous nucleation of one or several critical droplets which always appear at 
the system's border. That is, the free border acts as a droplet condenser. This is so be- 
cause it is energetically favorable for the droplet to nucleate at the border. Apart from 
the observed heterogeneous nucleation, the properties of the metastable- stable transition 
in equilibrium ferromagnetic nanoparticles do not change qualitatively as compared to 
the periodic boundary conditions case. [7] In our nonequilibrium system we observe a 
similar behavior, namely heterogeneous nucleation and the same qualitative nucleation 
properties. However, the fluctuations or noise that the nonequilibrium metastable sys- 
tem shows as it evolves towards the stable state subject to the combined action of free 
borders and the nonequilibrium perturbation are quite unexpected. 



As illustrated in Fig. 3. a, the relaxation of magnetization occurs via a sequence of 
well-defined abrupt jumps. That is, when the system relaxation is observed after each 
MCSS, which corresponds to a 'macroscopic' time scale, stricktly monotonic changes 
of m{t) can be identified that we shall call avalanches in the following. To be precise, 
consider the avalanche beginning at time t a , when the system magnetization is m(t a ), 
and finishing at %■ We define its size and lifetime or duration, respectively, as A m = 
\m(tb) — m(t a )\ and A t = \t{, — t a \. Our interest is on the histograms P(A m ) and P(A t ). 
Fig. 3.b shows the large avalanche size distribution P(A m ) for particle sizes R = 30 
and R = 60, after the extrinsic noise[8] (i.e. the trivial, exponentially distributed small 
avalanches) has been substracted. A power law behavior, followed by a cutoff is clearly 

observed. The measured power law exponents, P{A m ) ~ A m T< " R \ show size-dependent 
corrections to scaling. Similar corrections have been also found in real experimental 
systems.[9] In particular we find x(R = 30) = 2.76(2) and x(R = 60) = 2.06(2). The 

lifetime distribution also shows power law behavior, P(A t ) ~ A, y ' with a cutoff. 
Here we measure a(R = 30) = 3.70(2) and a(R = 60) = 2.85(2). This power-law 
behavior implies that avalanches are scale-free (up to certain maximum size and lifetime) 
in our nonequilibrium ferromagnet subjected to open boundary conditions. We also 
measured avalanches for p = in the circular magnetic particle case, and for p ^ 
in the periodic boundary conditions system. In both cases only small avalanches occur 
and the distributions are exponential-like, thus indicating the absence of scale invariance. 
That is, the combined action of free boundaries and impurities is behind the large, scale- 
free avalanches and essentially differs from the standard bulk noise driving the system 
and causing small, exponentially distributed avalanches only. The physical origin of this 
scale invariant behavior will be studied in a forthcoming paper. 

Summing up, in this paper we present preliminar results on the metastable behavior of 
a nonequilibrium ferromagnetic system. The presence of nonequilibrium conditions con- 
siderably enriches the observed phenomenology. In particular, we study the metastable- 
state mean lifetime for a lattice with periodic boundary conditions. Under the action of 
the nonequilibrium perturbation parametrized by p, the lifetime x{T,p, h) shows a max- 
imum as a function of T for certain nonzero temperature T max (p), then decreasing for 
lower temperatures. This counter- intuitive behavior, not observed in equilibrium, has 
some practical implications for real magnetic systems with impure ferromagnetic do- 
mains. We also observe that, under the action of both the nonequilibrium impurity and 
free borders, the metastable- stable transition proceeds by avalanches. These are power- 
law distributed, thus showing scale invariance (up to certain cutoffs). The chances are 
that our observations about the effect of the nonequilibrium conditions on the proper- 
ties of metastable states and their decay, which we can prove in our model cases, are a 
general feature of similar phenomena in real magnetic domains. 
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